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ABSTRACT 

I examine the role of dust grains in determining the structure of steady, cold, oblique 
C-type shocks in dense molecular gas. Gas pressure, the inertia of the charged compo- 
nents, and changes in ionization are neglected. The grain charge and rate coefficients 
for electron-neutral and grain-neutral elastic scattering are assumed constant at val- 
ues appropriate to the shock interior. An MRN size distribution is accounted for by 
estimating an effective grain abundance and Hall parameter for single-size grains. 

A one-parameter family of intermediate shocks exists for each shock speed v s 
between the intermediate signal speed vacos9 and y/2vACot9, where va is the pre- 
shock Alfven speed and 9 is the angle between the pre-shock magnetic field and the 
normal to the shock front. In addition, there is a unique fast shock for each v s > va- 

If the pre-shock density ?t-h ?Z> 10 5 cm -3 and the pre-shock magnetic field satisfies 
i?(mG)/nn(10 5 cm" 3 ) ^ 1 grains are partially decoupled from the magnetic field 
and the field and velocity components within fast shocks do not lie in the plane 
containing the pre-shock field and the shock normal. The resulting shock structure is 
significantly thinner than in models that do not take this into account. Existing models 
systematically underestimate the grain-neutral drift speed and the heating rate within 
the shock front. At densities in excess of 10 8 cm -3 these effects may be reduced by the 
nearly-equal abundances of positive and negative grains. 

Key words: magnetohydrodynamics - shock waves - dust grains. 
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1 INTRODUCTION 

The structure of shock waves in molecular clouds is deter- 
mined by the coupling between the magnetic field and the 
weakly ionized pre-shock gas. For shock speeds below 40-50 
kms" 1 (McKee, Chernoff & Hollenbach 1984) Lorentz forces 
within the shock front push the charged species through the 
neutrals, and the resulting collisions accelerate, compress, 
and heat the ambient gas. This process is slow because 
the charged species are rare, so the gas is able to radiate 
away a significant fraction of the heat while still within the 
shock front (Mullan 1971; Draine 1980), primarily through 
the emission of radiation in molecular rotational and vibra- 
tional transitions. Shocks of this nature are denoted 'C-type' 
(Draine 1980), or 'C*-type' if the gas is decelerated through 
a sonic point in the reference frame comoving with the shock 
front (Chernoff 1987, Roberge & Draine 1990). At higher 
shock speeds, the molecular coolants are dissociated, cooling 
can no longer keep pace with heating, and the gas pressure 
becomes dynamically significant. A thin, viscous sub-shock 
forms within the front, and the shock is termed 'J-type' 
(Draine 1980). 

C-type shocks efficiently convert the heat produced by 



collisions within the shock front into molecular line emis- 
sion, and are responsible for much of the intense infrared 
H2 and CO line emission observed towards the Orion-KL 
region (Draine & Roberge 1982; Chernoff, McKee & Hol- 
lenbach 1982; Smith & Brand 1990; Smith, Brand & Moor- 
house 1991; Chrysostomou et al. 1997). The warm molecu- 
lar environment within C-shocks can drive molecular chem- 
istry such as the conversion of atomic oxygen into OH 
and then into water (Draine, Roberge & Dalgarno 1983; 
Kaufman & Neufeld 1996a, b). In addition, the large drift 
speeds can drive exothermic ion-neutral chemistry within 
the shock front (Flower, Pineau des Forets & Hartquist 1985; 
Pineau des Forets, Flower & Hartquist 1986; Draine & Katz 
1986a,b). 

Shock models have generally been forced to be 'copla- 
nar', that is, the drift velocities and magnetic field within the 
shock front are assumed to lie in the 'shock plane', the plane 
containing the pre-shock magnetic field and the shock nor- 
mal (see Fig. 1). This holds for fastj^] shocks if the charged 



* 'Fast' and 'slow' have occasionally been used in the astrophysi- 
cal literature to denote J-type and C-type shocks respectively. In 
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species are tied to the field lines, but charged grains drift 
obliquely to the shock plane if they become partially de- 
coupled from the field lines by collisions with the neutrals 
(Draine 1980). However, for the sake of simplicity, the com- 
ponent of the grain drag force on the neutrals perpendicular 
to the shock plane has usually been supressed (e.g. Draine et 
al. 1983; Wardle & Draine 1987). This approximation breaks 
down for high pre-shock densities, partly because the frac- 
tional ionization of the gas decreases, permitting the grains 
to make a larger fractional contribution to the drag on the 
neutral gas, but also because the grains become poorly cou- 
pled to the magnetic field and the grain drag vector is further 
tilted out of the shock plane. 

This motivated calculations of the structure of grain- 
dominated C-shocks in which the vector components per- 
pendicular to the shock plane were retained. Pilipp, 
Hartquist & Havnes (1990) showed that the thickness of 
perpendicular shocks was significantly changed for pre-shock 
densities in excess of 10 7 cm -3 . In this case, the symmetries 
of the geometry constrain the magnetic field direction to be 
the same throughout the shock front, even though the drift 
velocities need not lie within the shock plane. 

Oblique shocks introduce an important extra degree of 
freedom, in that B±, the magnetic field component per- 
pendicular to the shock normal, may rotate within the 
shock front before returning to the shock plane downstream. 
Oblique shocks were studied by Pilipp & Hartquist (1994), 
who found solutions for shock speeds below a few kms - 
in which B± rotates by 180° (in either sense) within the 
shock front. Pilipp & Hartquist were unable to find solu- 
tions for shock velocities in excess of a few kms -1 , and 
suggested that steady shock solutions may not exist in this 
regime. This is disconcerting, as shocks with speeds in excess 
of 20 kms -1 are required to model the observed emission 
from a variety of sources. In addition, the development of a 
sound theoretical understanding of magnetohydrodynamics 
of dense, dusty gas is important for scenarios for the forma- 
tion and collapse of molecular cloud cores and the formation 
of protostellar discs (e.g. Nishi, Nakano & Umebayashi 1991; 
Ciolek & Mouschovias 1993; Neufeld & Hollenbach 1994; Li 
& McKee 1996; see also the review by Hartquist, Pilipp & 
Havnes 1997). 

The behaviour of B± in the solutions obtained by Pilipp 
& Hartquist is fundamentally different to the earlier copla- 
nar oblique shock models (Wardle & Draine 1987; War- 
dle 1991a; Smith 1992), in which the pre-shock and post- 
shock B± vectors are parallel rather than antiparallel. The 
upstream and downstream states across magnetohydrody- 
namic (MHD) shocks are related by jump conditions that 
are independent of the detailed nature and behaviour of the 
fluid, and therefore of the actual shock structure. In partic- 
ular, the transverse magnetic field increases, changes sign, 
or decreases across the fast, intermediate, or slow shocks re- 
spectively (e.g. Cowling 1976; Kennel, Blandford & Coppi 
1989). This implies that the solutions found by Pilipp & 
Hartquist are C-type intermediate shocks, whereas previ- 
ously studied C-shock solutions are fast shocks. The jump 
conditions show that the downstream state for intermedi- 



this paper I use 'fast' and 'slow' in the traditional MHD sense, 
referring to the fast and slow signal speeds. 



ate shocks is unphysical for shock speeds a little above the 
Alfven speed, and this explains the breakdown of the solu- 
tions for mildly super- Alfvenic shocks. Further, for a given 
shock speed, Pilipp & Hartquist found a one-parameter fam- 
ily of shock solutions rather than the single shock solution 
found in previous studies, and noted that this has also been 
found to be the case for intermediate shocks in resistive 
MHD (Wu 1988a; Kennel, Blandford & Wu 1990). 

Having identified these solutions as intermediate shocks, 
there still remains the important issue of the existence of fast 
shock solutions. On grounds of physical continuity, one ex- 
pects that adding grains to a coplanar, fast C-type shock 
should distort the structure out of the shock plane rather 
than do away with it entirely. Indeed, there are hints that 
steady shock solutions may exist. Following standard prac- 
tice in C-shock modelling, Pilipp & Hartquist derived a set of 
ODE's describing the shock structure, and integrated them 
from a perturbed upstream state towards the shock front. 
There is, however, a difference in the number of free param- 
eters describing the initial perturbed state. In the coplanar 
case there is one free parameter - the amplitude of the ini- 
tial perturbation in the pre-shock fluid (Draine 1980; War- 
dle & Draine 1987). Physically, this specifies the point in 
the shock precursor at which the integration towards the 
shock front begins, so there is a single shock solution for a 
given shock speed, the fast shock. Dropping the assumption 
of coplanarity introduces a second parameter, which can be 
regarded as the degree of rotation of the perturbation out of 
the shock plane at the initial point. The freedom introduced 
by the new parameter potentially allows a family of solutions 
to exist (the intermediate shocks), and Pilipp & Hartquist 
found that there is a critical value at which the sense of rota- 
tion of B± across the intermediate shock solutions changes 
sign. This suggests that the fast shock solution corresponds 
to this critical value, but that the solution cannot be found 
by integration from the perturbed upstream state because 
of finite numerical precision. 

This paper demonstrates that fast non-coplanar C- 
shock do exist. As the issues addressed here are fundamen- 
tal in nature, I keep the formulation of the multifluid shock 
problem, described in the next section, as simple as possi- 
ble. In particular, I assume that the molecular gas consists 
of four distinct species: neutrals, positive ions, electrons or 
PAHs, and negatively charged grains. Gas pressure is ne- 
glected, and ionization balance and chemistry are also ig- 
nored. These simplifications allow all of the physical quan- 
tities in the shock to be expressed in terms of the com- 
ponents, B x and B y , of B±. Thus a shock solution can be 
represented by a plot of B y (z) vs. B x (z), where z is the coor- 
dinate along the shock normal. This two-dimensional phase 
space, explored in is a powerful tool for examining the set 
of shock solutions given the pre-shock conditions and shock 
speed. In §^ I show that both fast and intermediate shock 
solutions exist for low shock speeds, and that the topology 
of the phase space prevents the fast shock solution from be- 
ing found by integration from upstream to downstream. At 
higher speeds the intermediate solutions become unphysical, 
but the fast solutions persist. I show that marginally cou- 
pled grains dominate the grain drag, and that their effect on 
the shock structure is dramatic, and examine the effects of 
an MRN grain-size distribution. The implications of these 
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results are discussed in and a summary is presented in 



2 FORMULATION 

The fluid is assumed to be weakly ionized, in the sense 
that the inertia and thermal pressure of the charged compo- 
nents, and the collisional coupling between charged species, 
are unimportant. This is generally an excellent approxima- 
tion in molecular clouds (Shu 1983; Mouschovias 1987), and 
within C-type shock waves as the fractional ionization in 
dense clouds is less than 1 part in 10 6 . Processes that create 
and destroy different species within the shock front, i.e. ion- 
ization, recombination, and chemistry, are important (Pilipp 
et al. 1990), but will be neglected here for simplicity (see gq). 
Further, the pressure in the neutral component is set to zero, 
a justifiable approximation in C-type shocks where the dy- 
namics of the flow is weakly dependent on the gas pressure 
because the neutral component is supersonic throughout the 
shock front (Wardle 1991b). 

Although the fluid temperatures do not directly influ- 
ence the dynamics, they affect the rate coefficients for elastic 
scattering between the charged species and the neutrals (see 
IpH). In addition, the electron temperature affects the cou- 
pling of the grains to the magnetic field because the grain 
charge is determined by the sticking of electrons to grain 
surfaces (eq. ^5|). However, I assume that the rate coef- 
ficients and the grain charge are constant throughout the 
shock front (with values chosen to be representative of the 
shock interior) so that equations for the temperatures need 
not be integrated. 

The equations describing shock structure developed in 
are therefore generally simpler than those of Pilipp & 



2.1 



Hartquist (1994), who explicitly followed the fluid tempera- 
tures, grain charging, and chemistry. They are more general 
in allowing an arbitrary number of charged species, and in 
not assuming that the magnetic flux is frozen into the elec- 
tron fluid. The latter has little effect for the densities of order 
10 6 cm -3 relevant here, but will allow the equations to be 
applied at higher densities (in excess of 10 9 cm -3 ) that may 
be relevant for models of shocks associated with maser re- 
gions or of accretion shocks onto protostellar discs (Neufeld 
& Hollenbach 1994) where the electron density is low, and 
ions and PAHs become decoupled from the magnetic field. 

The adopted properties of the charged species present 
in the gas - ions, electro ns, negatively charged grains and 
PAHs - are discussed in E 2.2 . 



2.1 The equations for shock structure 

Each component of the fluid is characterised by particle 
mass m and charge Ze, mass and number densities p and n, 
and velocity v. Unsubscripted quantities refer to the neutral 
fluid, and subscripts i, e and g refer to the ions, electrons, 
and negatively charged grains respectively. The generic sub- 
script l j' will be used to denote any charged species, i.e. 
j G {i,e,g}, and the equations are trivially generalised to 
include additional charged species by expanding this set. 

The shock front is assumed to be steady and plane- 
parallel, propagating at speed v s into a uniform medium of 
neutral density po permeated by a magnetic field Bo which 
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Figure 1. A sketch of the Cartesian reference frame used in this 
paper. In this frame the shock is at rest, and the unshocked ma- 
terial moves into the shock along the z-axis at speed v s . The 
upstream magnetic field makes an angle 9 to the incoming veloc- 
ity, and lies in the x-z plane, which is referred to in the text as 
the 'shock plane'. 



makes an angle to the shock normal. We adopt a Cartesian 
coordinate system (x, y, z) in the shock frame (see Figure 1) 
with the z-axis parallel to the direction of shock propagation, 
and the pre-shock magnetic field lying in the x-z plane. The 
incoming pre-shock gas then has velocity v s z. The continuity 
equations are 



pV z = pOVs 

for the neutral fluid, and 



(1) 



PjV jz = PjQVs (2) 

for each of the charged species. The momentum equation for 
the neutral component of the fluid is: 

d 



dz 



{pv z v) 



E 

3 



(3) 



where 7j is related to the rate coefficient for elastic scattering 
between particles of species j and the neutrals, < av >j, by 

< (TV >j 



7i 



m, + m 



(4) 



The drift of each charged species through the neutral gas 
is determined by a balance between electromagnetic forces 
and collisional drag (e.g. Shu 1983): 



,Z s e (e + x fi) + j jPjP (v -Vj)=0. 



(5) 



The inertia of the charged species has been neglected, which 
is valid provided that the time scale over which conditions 
change (i.e. the flow time through the shock front) is long 
compared to the collisional time with the neutrals. Equation 
(fl) can be rewritten as 



, a (B ■ E')B 



0$ E' X B 



B 3 ' 1 + /3] 



B 2 



+ 



Pj B x (E' x B) 

1 + P] B 3 



(6) 



where 

E' = E + vxB/c (7) 

is the electric field in the frame comoving with the neutrals, 
and Pj is the Hall parameter for species j: 

Zi&B 1 



P 



rrijC ~jjp 



(«) 



4 AT. W-W/e 



i.e. the product of the gyrofrequency and the time scale for 
momentum exchange with the neutral fluid. Note that the 
sign of j3 j is determined by Zj . The magnitude of the Hall pa- 
rameters determines the magnetic flux transport in a weakly 
ionized medium (Draine 1980; Konigl 1989). When \j3j\ > 1, 
the electric and magnetic stresses individually dominate the 
drag force in equation (5) and must almost cancel, thus 
E ss —(vj x B)/c and the particles are tied to the mag- 
netic field lines. In the opposite limit, collisions with the 
neutrals decouple the charged particles from the magnetic 
field lines and Vj ~ v. 

Rather than explicitly consider the drift velocity of each 
charged species, it is often more convenient to consider the 
expression for the current density 



J = e^rijZjVj 



0) 



obtained by using (|6j) to substitute for Vj\ 

J = a {B ■ E')B + o-iB X E' + a 2 (B X E') X B , (10) 

where the components of the conductivity tensor are the 
conductivity parallel to the magnetic field, 



ec v— * 

j 

the Hall conductivity, 
ec \ - rijZj 



o\ = 



B £—i 1 + /?? 



and the Pedersen conductivity 
ec \ - UjZjPj 



a 2 



b z^ i + a? 



(ii) 



(12) 



(13) 



(Cowling 1957; Nakano & Umebayashi 1986). 

Now consider Maxwell's equations. The induction equa- 
tion yields VxE = 0, implying that E x and E y are constant. 
Evaluating them far upstream gives: 



E X = Q, 
and 

E y = Box ■ 

c 

Finally, V • B — implies that 
B z = Bq z . 



(14) 
(15) 
(16) 



Gauss' law for the electric field is replaced by the condition 
for charge neutrality: 



Ampere's law reduces to 

dB x _ 47T 

dz c 
dB y _ 47T 

dz c 

and 



Jz = 0. 



(17) 

(18) 
(19) 

(20) 



Equation (^0|) is not independent as eqs (Q) and (|lj) guar- 
antee J z = upon noting that Vj z = v s far upstream. 

The derivatives in equation (^|) can be eliminated us- 
ing (j|), @, © and ©-© to obtain expressions for 
the conservation of the x, y, and z components of the total 
momentum flux: 



POVsVy 

and 



B x Bz 
4n 



4tt 



B() x Bz 

4-k 



(21) 
(22) 



PoVsVz + 



B 



2 + Bl 



= POVs + 



(23) 



Bq x 

8tt r '"" s ' 8tt ' 

Thus we have a s yste m of two ordinary differential equa- 
tions in B x and By, (JlSJ) and ( P.9L along with the algebraic 
relations ©,(§), and 

2.2 Charged species 

The pre-shock abundances of charged species are determined 
by a balance between cosmic-ray ionization and recombina- 
tion occuring in the gas phase or on grain surfaces, and are 
therefore sensitive to the adopted grain model (Nishi et al. 
1991). 

When grains are of a single size a J> 0.1 p, recombina- 
tion occurs primarily in the gas phase for tih ^ 10 10 cm' 3 , 
yielding a fractional ionization 



i*. « io- 8 f 



mi 



10 6 < 



-1/2 



(24) 



and an almost identical electron abundance (Elmegreen 
1979; Nakano & Umebayashi 1986). The grains in the pre- 
shock gas carry a single negative charge, and assuming the 
the total mass of grains is 0.01 that of hydrogen, the grain 
abundance is n g /riH = 2.5 x 10~ for a — 0A/i (Draine et al. 
1983), or 1.6 x 10 -12 for a = O.lp. Under these conditions, 
the charging of the grains within the shock front proceeds by 
the sticking of electrons from the gas phase, and is limited by 
the Coulomb repulsion of incoming electrons by the grain. 
The resulting mean charge depends on the grain radius and 
the electron temperature (Draine 1980): 

- ikT " -(^k)(^. «> 



-4.8 



valid for \Z g \ 3> 1. The charges of individual grains fluctuate 
about the mean, but at any instant most are within 1 unit 
of the mean (Elmegreen 1979). 

However, grains possess a power-law grain-size distribu- 
tion (Mathis, Rumpl & Nordsieck 1978; Draine & Lee 1984), 

dn A -3-6 /or\ 
— — = Anna (26) 

da 

where n^a) is the number density of grains with radii smaller 
than a, A — 1.5 x 10~ 25 cm~ 2,5 , and the distribution runs 
from grain radii between lower and upper bounds ai = 50A 
and 02 = 2500A. The lower bound is poorly constrained 
by observations, and the upper bound is uncertain in dense 
clouds, where grains are enlarged by icy mantles. In addition 
PAHs with a ~ 3A and n/n H « 2 x 10~ 7 (Leger & Puget 
1984; Puget & Leger 1989), may be present, either as a dis- 
tinct population or as the tail of an MRN distribution with 
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a\ m 3A. For the sake of definiteness, I shall treat the PAHs 
as a distinct population. The MRN size distribution and the 
presence of PAHs increase the grain surface area and reduces 
the fractional ionization in the gas and the electron density 
by several orders of magnitude (Nishi et al. 1991; Kaufman 
& Neufeld 1996a). The electron density is comparable with, 
or less than, the density of grains and the charging of grains 
by electron sticking within the shock front is inhibited. Thus 
the grain charge remains at Z g = — 1, independent of T e in 
this case. 

These two cases are considered here. First, for compari- 
son with Pilipp & Hartquist (1994), I consider large grains of 
a single size. The charged species in this case are ions, elec- 
trons and negatively charged grains. The effect of an MRN 
grain-size distribution and PAHs are also considered. In this 
case the charged species are ions, electrons or PAHs, and 
negatively charged grains. The grains are still considered to 
be of a single size, but the grain size and abundance are cho- 
sen so that the contribution to the conductivity of the gas 
appr oximates that of the entire grain-size distribution (see 
& 

I assume that the ions are singly charged and have mass 
30mH, as is typical for molecular or metal ions (e.g. HCO + , 
Mg + ). Grains are taken to be spherical, all of the same ra- 
dius a, with internal density 2.5 gem -3 and total mass 0.01 
of the mass in hydrogen. 

The rate coefficient for elastic scattering with the neu- 
trals, < ov > j , depends in general on the fluid temperatures 
and the drift speed \vj — v\. The scattering cross-section for 
ions has a l/« dependence for ions with drift speeds below 
20 kms -1 , and the rate coefficient for ion-neutral scattering 
is constant, < av >i~ 1.6 x 10 _9 cm 3 s _1 , yielding a Hall 
parameter 



Pi « 4600 



(lmc) ( 



10 6 < 



(27) 



The PAHs have geometric cross-sections na ~ 3 x 
10 _16 (a/A) 2 cm -2 , of order the polarization cross-section for 
elastic scattering from neutrals, 1.6x 10~ lj x 10kms _1 /|t;g — 
v\. The PAHs are therefore treated as negatively charged 
ions, with < av >pah = < crv >;, and /?pah = — Pi- 

The elastic scattering cross-section for grains is the ge- 
ometric cross-section, taken to be a g = tyo 2 , while for elec- 
trons, a e ~ 1 x 10~ 15 cm 2 at energies below a few eV (Gi- 
lardini 1972). As Oj is energy- independent for both these 
species, the rate coefficient may be written as <JjUj where 
the effective velocity is 



128 
~9n 



kT kTj 



+ \v.i - v 



1/2 



(28) 



(Draine 1986). The rate coefficient for electron scattering is 
determined by the electron thermal velocity, which is far in 
excess of both the electron-neutral drift speed and the neu- 
tral thermal velocity. The effective relative velocity between 
electrons and neutrals may therefore be written 



/ 128fcT e \ 1/2 



V 9nm e 

and the electron Hall parameter is 
6 / B \ ( n H 



(29) 



Pe « -2.5 X 10' 



(lmc) (i 



,10 6 cm- 3 
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The grain thermal velocity is negligible compared with that 
for the neutrals, and 



/ 128fcT 



+ \ v s - «l 



1/2 



(31) 



V 9nm 

The grain Hall parameter is then [eq. (35) of Draine (1980) 
B 



io- 



The rate coefficient for grain-neutral scattering depends 
on the neutral gas temperature and the grain-neutral drift 
speed. The grain charge and the electron scattering co- 
efficients depend on the electron temperature. A calcula- 
tion of the fluid temperatures is beyond the scope of the 
present work, so we assume that u e , u g and Z g are con- 
stant throughout the flow, being fixed at values representa- 
tive of conditions within the shock front. In particular, for 
the weak shock models (v 3 /va ~ 1.5) I adopt T = 100 K, 
T e — 200 K, and \v — v g \ — lkms -1 ; and for the strong 
shock models (v s /v A = 10), T = 1000 K, T e = 2000 K, and 
\v — v„ = lOkms -1 . The grain charge is either calculated 
from (I25J) using the adopted values of T e or is set to — 1, 
depending on the chosen grain model. 



3 ANALYSIS 

The shock structure is determined by the pair of differen- 
tial equations ( |l8[ ) and ( |l9| ) for B x and B y as a function of 
z, and the associated algebraic equations. There is no ex- 
plicit reference to z, so given a choice of parameters, each 
shock solution can be represented as an integral curve in 
a phase space with coordinates B x and B y , i.e. as a plot of 
B y (z) vs B x (z). This proves to be extremely useful in under- 
standing the relationship between different shock solutions. 
Trajectories representing shock solutions begin and termi- 
nate at points representing upstream or downstream states, 
for which 

d 



d_ 

dz 



(By) = 



(33) 



dz 

( i.e. stationary points). The jump conditions, discussed in 



show that there are three stationary points, all lying 
on the B x axis, which are denoted by U, F, and I in increas- 
ing order of \B X \ (see Fig. Q). These points represent the up- 
stream state, and possible downstream states corresponding 
to the fast and intermediate shocks respectively. The topol- 
ogy of the trajectories near t he st ationary points is obtained 
through a linear analysis in § |3.2| , where it is shown that U, 
F, and I are a source, saddle, and sink respectively. Thus it 
is not possible to find a fast shock trajectory by starting at 
the upstream state and integrating through the shock front 
towards the downstream state F; the trajectory should be 
integrated in the reverse direction from F to U. 

To actually integrate equations ([l8]) and (jl9|) , the x and 
y components of the current density J must be determined 
from the algebraic relations given B x and B y at some point 
in the shock front. The expressions for total momentum flux 
(pl|)-(p3|) allow v to be expressed in terms of B, the neutral 
continuity equation then yields p, and the Hall parameters 
Pj can then be calculated for each charged species. If E' z 
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were known then the charged particle drifts, and hence the 
current density, could be calculated from eq. (^). Actually 
E' z is determined from the charge neutrality condition (|l7j), 
by substituting the z-component of Q and using (g) to re- 
late pj and Vj Z (§3.3). The presence of N charged species 
allows TV — 1 solutions for E' z that yield drift velocities con- 
sistent with the charge neutrality condition at each point 
within the shock front. The choice of solution is dictated 
by requiring that the electric field vanish in the neutral rest 
frame ahead of (or behind) the shock front and be every- 
where continuous. 

Several regions of the phase space can be excluded on 
the grounds that the density of the neutrals, or of a charged 
species, becomes negative (§ 5.4). In particular, when B 2 be- 
comes too large, the ram pressure must become negative 
(see eq [p3[). The stationary point I becomes unphysical in 
this manner when the sh ock velocity is too large to allow 
intermediate shocks (see §3.1). When this is the case, there 
is also a locus in phase space where the ion and electron 
densities simultaneously diverge. However, it will turn out 
that the fast shock solutions avoid this locus. 

The effe ct of an MRN grain-size distribution is exam- 
ined in 



3.5 



where it is shown that it can be approximated 
by a single-size grain model with appropriately chosen 'ef- 
fective' grain Hall parameters and abundances. 

One interesting feature of the shock solutions is that 
within the shock front, the fluid motions and magnetic field 
are not restricted to lie in the plane containing the shock 
normal and the pre-shock magnetic field. The sense of mo- 
tion out of the plane (i.e. 'up' or 'down') is determined by the 
asymmetry in the elastic scattering properties of species of 
different sign (i.e.|ci| 7^ 0, see eq ( |12| ) ). Even if the Hall pa- 
rameters of all charged species satisfy \/3j\ ^ 1, fast shocks 
remain coplanar if the properties of negatively and posi- 
tively charged species are identical apart from the difference 
in sign. A cr iterion for this effect to be important is devel- 
oped in §3.6. In essence, there should be 'sufficient' poorly- 
coupled charged particles, with the charge on the particles 
being preferentially of one sign. 

Finally, in §3.6 I discuss the dimensionless parameters 
that determine the shock structure, and outline how the 
calculations presented in §W were carried out. 



3.1 Jump conditions 

The boundary conditions for equations ([ij|) and (|l§) are 
that the derivatives of B x and B y should vanish far away 
from the shock front. This implies that the current density, 
and hence (via eq Jl^]) the electric field in the rest frame 
of the neutrals also vanish at z = ±00. In turn, Q then 
guarantees that Vj = v far upstream and downstream of 
the shock front. This is, of course, expected, and shows that 
the MHD jump conditions which relate the upstream and 
downstream states can be obtained by regarding the weakly 
ionized fluid as a plasma of density p and velocity v, and 
equating the fluxes of mass and momentum of the combined 
fluid on either side of the front. The ionized component does 
not contribute significantly to the fluid density or mass flux, 
but does contribute to the momentum flux through the mag- 
netic field (see eqs |^l[-[^3|). An energy condition is not 
required, as it is assumed that T = 0. 

In general, the MHD jump conditions imply that the 



pre-shock and post-shock magnetic fields and velocities are 
coplanar, and three broad classes of shocks can be distin- 
guished on the basis of the behaviour of the field component 
transverse to the shock velocity, which either increases, de- 
creases, or changes sign in going from upstream to down- 
stream in fast, slow, and intermediate shocks respectively 
(e.g. Cowling 1976; Kennel et al. 1989). 

The signal speeds in the fluid play a key role in deter- 
mining the different kinds of shock the medium can sup- 
port. The combined fluid supports fast, intermediate and 
slow waves with signal speeds for propagation in the direc- 
tion normal to the shock front of 



B 2 \ 1/2 
v **=\4^) = VA 

where va is the Alfven speed, 

( B *Y /2 

Vint = ~ = VA COS 6 , 



and 



Vslo 



0. 



(34) 



(35) 



(36) 



respectively. The assumption that the plasma is cold elimi- 
nates shock transitions that have downstream velocities be- 
low the slow speed (i.e. slow shocks and three classes of 
intermediate shock, see Kennel et al. 1989). 
Far away from the shock front, 



E = -v/c X B 
which implies that 

VyB Z — V Z By = , 

and 

v z B x — v x B z = VsBq 



(37) 



(38) 



(39) 

Equations (^) and (pl[)-(|23|) for mass and momentum flux 
conservation also hold and combined with ( p8[ ) and ( pis| ) yield 
a complete set of jump conditions. 

Equations (^2|) and (53) both express v y downstream in 
terms of B y , implying that either B y — 0, or that 



4npov s 



(40) 



In the latter case, the z-components of the pre-shock and 
post-shock velocities are equal to the intermediate speed, 
and the transverse components of the magnetic field and 
fluid velocity change direction in such a way as to conserve 
the transverse momentum flux. There is no compression 
or entropy increase across such a transition, which exists 
in ideal MHD as a 'rotational discontinuity' (e.g. Cowling 
1976). In the non-ideal case it is unsteady, broadening dif- 
fusively with time (Wu 1988b), and we shall not discuss it 
further here. 

Thus we set B y = and v y = and use ( E l |) and ( |23[ ) 
to eliminate v x and v z in favour of B x in (fcS9])7 to obtain 
a cubic in B x whose solutions satisfy the jump conditions. 
The upstream state is a solution, so (B x — Bq x ) factors out 
of the cubic to yield a quadratic 



B 2 + B 0x B x + 2(B 2 - 4np i 







(41) 
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V s /\cosB 

Figure 2. A polar plot of the conditions for the existence of inter- 
mediate shocks. A point on the plane represents a shock velocity, 
with polar radius representing the shock speed v s , and polar an- 
gle representing the angle 8 between the pre-shock magnetic field 
and the shock normal. Intermediate shocks do not exist in the 
shaded region. 



with positive and negative roots that are candidate down- 
stream states for fast and intermediate shock transitions re- 
spectively. 

The roots correspond to physically acceptable down- 
stream states if there is compression across the shock front, 
i.e. for < v z <v s . From eq (E3|), the condition on B x is: 



1 < 



\BoxJ 



2 , , 8typqv 2 s 
R 2 



(42) 



When v s /va > 1, the positive root yields a fast shock, and 
when cos 8 < v s /va < V% cot 9 the negative root yields an 
intermediate shock (c.f. Kennel et al. 1989). The regime in 
which intermediate shocks are allowed is sketched in Fig. 2. 



3.2 Classification of the stationary points 

The stationary points can be classified as source, sink or sad- 
dle by linearizing the equations for shock structure around 
each point. I use eqs. (1) and (16)-(18) for the conservation 
of neutral mass flux and total momentum flux, but rather 
than explicitly considering the drift velocities of the charged 
species given by (^), it is more convenient to use (|lo|). Since 
E' — outside the shock front, the linearised version of this 
equation is 

SJ = a {B ■ 5E')B + diB X SE' + a 2 (B X 5E') X B (43) 

Assuming that the perturbations are steady, independent of 
x and y, and vary exponentially in z as e Az , the resulting 
dispersion relation is 



(lis — v\ cos 2 9)p ■ 



(cr 2 + cr 2 ) sin 2 9 + <too"2 cos 2 9 



[(u 2 — «a)A» — cr 2 ] + cr 2 cos 2 = 0, 
where 

4ttoa (cr? + erf) 



/'• 



c 2 Vs A 



(44) 



(45) 



The nature of the stationary point is determined by the sign 
of the real part of A, which is the same as that of fi. Noting 



that cro, cr 2 > 0, and using v a > va, we find that the real 
parts of the roots of the dispersion relation are positive, so 
the upstream state is a source. 

The dispersion relation can be applied to the down- 
stream states F and I by noting that a boost from the 
shock frame to a frame in which the shock is stationary 
and the downstream fluid velocity is normal to the shock 
does not affect the assumptions of steady state and the e Az 
dependence of the perturbations. Thus ( fl4| ) can be applied 
downstream on replacing v 3 by v z , va by the downstream 
Alfven speed, and 9 by the value of cos - (B x /B) far down- 
stream. Re-evaluating the dispersion relation at F (for which 
va > Vz > va cos 9) shows that the roots there are real and 
of opposite signs, and at I (where v z < vacosO) the roots 
have negative real part. Thus F is a saddle point, and I is a 
sink. 

Finally, we note that at U and I, when 



/ 2 



„2 1 '2 
2 . CTi + cr 2 

v A ) — -o— 



2 CT2 



+2(« 2 



2 \l/2/ 2 

■v A ) {Vz 



2 2 

va cos ( 



1/2 cos6> |ai 
sin 2 9 cro 



> (46) 



the roots of the dispersion relation, eq. ([m]) are complex con- 
jugates and trajectories spiral out of U and into I. Typically, 
ctq S> o\ + cr| and v 2 v\, and the criterion becomes 



cos 9 

CT 1 > ~ 2~a ^ 2 I 

sin 9 

or \/3 g \ £ 0.3 for 



(47) 



30° 



3.3 



The charge neutrality condition 

To actually integrate equations (^) and ([L9|), the x and 
y components of the current density J must be determined 
given B x and B y . This is achieved as follows. The expressions 
for total momentum flux (^])-(^3|) allow v to be expressed in 
terms of B, the continuity equation then yields p, and thus 
the Hall parameters j3j can be calculated for each charged 
species. At this point, E' x and E' y can be found from (jjj), (jlij) 
and dl5|), so that if E' z (or equivalently E z ) were known then 
the charged particle drifts, and hence the current density, 
could be calculated from eq. (p|). Actually E' z is determined 
by the charge neutrality condition (|T^), by substituting the 
z-component of (ph and using (□) to relate pj and Vj Z : 



n j0 Zj 



Noting that eq. (£ 
form 



(48) 



implies that Vj z can be written in the 



v jz = n j0 Z J {p i E' z + q j ) 
where 



Pj = 
and 

1i = 



1 (3,(1 + /3 2 B 2 Z /B 2 ) 
n j0 Z 3 (1 + 13]) 

1 [ c0] (E'XB), 

n j0 ZjV s [ Vz 1+/3] B 2 

c(3j (B-E')B Z 



(49) 



(50) 



+ 



1 + /3 2 



(51) 
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eq. ( 

E 



jj) can be recast as 
0. 



PjE' z + qj 



(52) 



In this expression, pj and qj are known quantities, and the 
task is to determine E' z . If N is the number of charged 
species, the sum regarded as a function of E' z has N poles 
at E' z = Ej where Ej = —qj/pj- If the poles are ordered so 
that E\ < E2 . ■ ■ < E n , then there is a zero between each 
pair Ej, Ej-\. Thus locally there are several values of E' z 
which will produce particle drifts with no net current. The 
correct choice is dictated by the requirement that E' z (z) be 
smooth and go to zero ahead of and behind the shock. 

Multiplying (^) by J^-fe-E^ + qj) yields a polynomial 
in E' z which for N charged species has N-l solutions. For the 
case N = 3 relevant to this paper, eq. (^) yields a quadratic 
in EL: 



a 2 E' z +a 1 E' z +a =0, 



where the required root is 

e: = 



-ai 



\/ a\ — 4aofl2 



2a 2 



(53) 



(54) 



3.4 Forbidden regions of phase space 

Several regions in the B x -B y phase space are forbidden in 
the sense that the density of one or more species formally 
becomes negative. 

Firstly, the magnetic field strength is limited by : 



B x + B y 
B 2 



< 



2vi 



+ sin 



(55) 



otherwise eq ( |23[ ) implies that v z would need to be negative 
to conserve the total z-momentum flux. As the magnetic 
field approaches the limiting value, v z tends to zero and the 
neutral density diverges. 

An additional restriction is imposed by requiring that 
Vj Z > for each charged species. The locus Vj z = cor- 
responds to diverging particle density (see eq. and 
the charge neutrality requirement forces the density of an 
oppositely charged species to also diverge there. Thus for 
the system of ions, electrons, and negatively charged grains 
considered here, there are two loci: Vi Z = v ez — 0, and 
Viz = v g z = 0. Which of these loci is relevant is intimately 
related to the choice of E' z as determined by the charge neu- 
trality condition. As the ions and electrons are better tied 
to the field lines than the grains, the choice of E' z shows that 
it is the first condition that is relevant. For \f3 e \ and (3i 3> 1, 
the locus Vi Z = v e z — is given by 



Bo 



+ 



Ey 

B 



(56) 



a circle centred on a point on the negative B x axis. The 
region enclosed within the locus is unphysical, as the elec- 
tron and ion densities become negative there. This locus 
exists only if the right-hand side of eq. ( [56| ) is positive, i.e. 
when v s > \/2?MCot 9. This is precisely the criterion for the 



non-existence of the intermediate shocks, so the unphysical 
region in the B x -B y plane is only relevant for fast shock 
waves. 

The constraint that Vj Z > arises from the continuity 
equations for the charged species and the requirement that 
the charged particle densities should remain positive. To 
some degree this is an artifact of the neglect of ionization and 
recombination processes ocurring within the shock front, re- 
combinations serving to limit the divergence of the charged 
particle densities. In principle, charged particles could drift 
backwards at some point within the shock, with a steady 
state being maintained by recombinations around the points 
where Vj Z = 0. In any case, the increase in charged parti- 
cle density increases the currents that can flow within the 
shock front and prevent the trajectories from crossing into 
the unphysical zone (see §4.2). 



3.5 The grain-size distribution 

The majority of studies (e.g. Draine et al. 1983; Pilipp & 
Hartquist 1994) assumed grain sizes of 0.4 fi and 0.1 p, re- 
spectively. However, the smaller grains in an MRN distribu- 
tion dominate the conductivity of the gas (Nakano & Ume- 
bayashi 1986) and determine the shock structure (Kaufman 
& Neufeld 1996a). If all grains have the same drift veloc- 
ity, the cross-section for collisions with neutrals scales with 
grain radius as a 2 , and the differential contribution to the 
drag by grains with radii in the range a to a + da scales 
as a~ 15 da oc da~ 5 . (The decrease at larger radii will be 
a little steeper as grains become more strongly tied to the 
neutral gas.) The dust contribution to the drag force on the 
neutral gas is therefore dominated by the smallest grains, 
but a broad range of sizes contributes significantly to the 
total drag force. 

can be applied in this 



In principle, the equations of §2.1 



case with each size grain treated as a different species, and 
the sums over j converted to integrals over grain size. Such 
a treatment is complicated by the dependence of \v g —v\ on 
grain size, and lies beyond the scope of the present paper. 
Instead, I develop an approximate treatment in which an 
effective grain size and number density are determined. 

The rationale for selecting the 'effective' pre-shock grain 
no. density n g and Hall parameter j3 g is that the grain con- 
tribution to the components of the conductivity tensor per- 
pendicular to B, i.e. o 1 and 02, should match that for the 
entire size distribution. It is not possible to simultaneously 
match the contribution to the conductivity parallel to the 
magnetic field, ao, but this is dominated in any case by the 
species that have large Hall parameters, i.e. ions, and elec- 
trons or PAHs (see fill]). Thus n g and (3 g must satisfy 



a2 dn g /da 



da = Nh 



and 

n B \/3 g \ _ f a2 \P g \dn g /da 

/ a 1 



1+/?' J a , 



da = NI 2 



where 



N = ^Anna' 2 5 w 3.4 x 10~ iu n H 



r * v 



(57) 



(58) 



(59) 



V50A/ 
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is the total number of grains with radii between ai and ai 
apart from a negligible correction factor of [1 — (ai/a2) 2 ' 5 ], 



h = \Pl 
and 

h = fft 



5/4 



fl P 



1/4 



l+P 



^dp, 



|9i 



02 



ft/ 4 
l + P 2 



■dp, 



(60) 



(61) 



where ft and ft are the absolute values of the Hall param- 
eters for grains of sizes a\ and ai respectively. The integrals 
over o have been converted to integrals over oc a~ 2 . Thus 

ft = -h/h (62) 

and 

n g = (1 + p 2 g )hN . (63) 

The integrals for I\ and I2 are solved analytically be- 
low, but it is instructive to derive approximate values before 
presenting the exact results. These can be obtained by set- 
ting 1 + ft in @ and (|l|) to 1 for P < 1 and to p 2 for 
0>1. This yields 



h 



and 



h 



4„-5/4 
3 Pi 



pi <C ft < 1 

ft « K ft 



(aa/ax) 3 ^- 2 1 « ft « ft 



fft ft < ft < 1 
5ft" 1 1 < ft 



(64) 



(65) 



The approximation is poor for 7i when ft <C 1 <C ft, as 
there is a significant contribution for near 1. Thus the 
expression given above includes a correction factor of 1/2. 
These expressions yield estimates of the effective grain pa- 
rameters: 



Pi 

and 



2 

Ipi^ 



ft < ft < 1 
ft < K ft 
3(a 2 /a 1 )- 3 '' 2 1< ft < ft . 



(66) 



ft « ft < 1 

ft < K ft 
15(a 2 /ai)- 3/2 K ft < ft . 



(67) 



These approximations show that when all of the grains 
are decoupled from the magnetic field (ft <C 1), the effec- 
tive Hall parameter is roughly half that of the most strongly 
coupled grains, and all of the grains contribute to the non- 
coplanarity. When all of the grains are coupled (ft ^> 1), 
ft w 3(a2/ai) 1,/2 ft , and the number of participating grains 
is of order one per cent. However, the most relevant case 
is the intermediate regime, in which the smallest grains 
are coupled to the field, whereas the largest grains are not 
(ft < K ft), as from (El), 



ft w 8.6 



B/mG 



72h/10 6 cm~ 



V5CL47 



VlOkms- 1 / V ; 



and ft w 2500/3i . Then the effective Hall parameter is of or- 
der unity, and a tenth of the grains or more can be regarded 
as participating. 




lQ glO Pi 



Figure 3. The effective single-size grain parameters correspond- 
ing to an MRN distribution between 50 and 2500 A, plotted as 
a function of the magnitude of the Hall parameter of the 50 A 
grains, ft. The solid curve shows the effective grain no. density 
relative to the total grain no. density, and the dashed curve shows 
the effective Hall parameter ft relative to ft. ft is typically in 
the range 0.1 to 1 in dense clouds. 



7i and I2 can be found analytically by making the sub- 
stitution x = /3 1//4 and using the identity 



■ dx — 



n 

Esin mOk tan -1 ( ) 
L V sin6l fc / 



k=l 
1 



cos mOk ln(l — 2x cos 9k + x 



where 



2k 



2n 



(69) 



(70) 



which is valid for m and n natural numbers and m < 2n 
(Gradshteyn & Ryzhik 1965, 2.146). Then 



Ii = 5ft 



h = 5ft 



5/4 



5/4 



hAiP\ /A ) 



h,A {Pi 

IiAP^)-IiAPT)+pT-pT 



(71) 



(72) 



The corresponding values of n g and ft are plotted in 
Figure ^. As ft ranges from 1 to 10, n g /N decreases from 
about 1 to 0.5, and |ft| increases from 0.6 to 3. Thus the 
effective Hall parameter typically lies between 0.1 and 2, 
and the effective number density is between 0.5 and 1 of the 
actual charged grain number density. 



3.6 Criteria for non-coplanarity of fast shocks 

The general criterion is simply j<7i| <; 02. A rough criterion, 
assuming that ions and electrons are well-coupled is that 
the contribution of poorly-coupled grains to the current is 
a significant fraction of the current that would flow within 
the shock in the absence of grains. The former can be es- 
timated by assuming that the grain-ion drift is 1/2 of the 
shock speed, the latter from the magnetic field gradient in 
a perpendicular C-type shock, yielding: 
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\n a e\Z g \v s ^ 
where 

-^shock 



c V2v s /vaB 

4n I/shock 



a/2 



"a 



71i < (TV >i 



(73) 



(74) 



is the characteristic shock thickness in the absence of grains 
(Wardle 1990). This criterion can be rearranged to yield: 

^Ap^l. (75) 

A more rigorous derivation begins by considering the 
grain contribution to the net drag force on the neutrals, or 
equivalently, the Lorentz force per unit volume: 



J X B _ (Ti(B X E') X B + a 2 B X E' 



(76) 



Assuming that n 9 < n e ~ n;, and that \/3 e \ j3i S> 1, then 
the electron contribution is dominated by the ions, and the 
ion contribution to eq. (|76[) is 



The grains 

( J x B \ 

V ec /, 



Pi 

contribute a force 



(77) 



Jsj^ (BXE')XB + B x £'(78) 



which dominates the ion drag provided that 



riiZi 



(79) 



aA+A 

The grain contribution is non-coplanar when Bg Si 1, thus 
combining this requirement with ( |7S| ) yields (j75|). 

The criterion for non-coplanarity, eq. (75), is the oppo- 
site to that given by Pilipp & Hartquist (1994, §3). Their 
sense of the criterion implies significant rotation in the limit 
of no grains, n g — > 0, yet fast shock waves in this limit are 
planar. This was based upon a linear analysis of the shock 
precursor, and derived for the growing mode that is rele- 
vant to intermediate shock precursors, rather than the fast 
shocks which the criterion of eq. (75) addresses. Even when 
the Hall current vanishes - in the pure ambipolar diffusion 
limit - rotation in intermediate shock precursors is enforced 
by the requirement imposed by the jump conditions that 
the sign of B x changes across the shock front. Fast shocks, 
however, will not exhibit rotation unless the Hall current is 
important. 



3.7 Numerical solution 

The equations are put in dimensionless form by using v 3 , 
Bo, v 3 Bo/c as units of velocity, magnetic field strength and 
electric field respectively. The parameters that appear in the 
equations are 

(i) A = v s /va, the Alfven number, which determines the 
strength of the shock, 

(ii) 8, the angle between the upstream magnetic field and 
the shock velocity, 

(iii) Pjo, the Hall parameters in the pre-shock gas (eq. 
||), which determine whether the field is frozen into each 
species {13 j i; 1), and 



(iv) "/jPjo/^^jPjo, which is a relative measure of the 
contribution of species j to the net drag force on the neutral 
component given equal drift velocities of the charged species. 

The overall normalisation of the drag forces appears in the 
expression for a characteristic length scale, 



HPjO ■ 



(80) 



which appears as the unit of length in the dimensionless 
equations for dB x /dz and dB y /dz. Although this is approx- 
imately equal to the shock thickness L s hock when the charged 
species are well coupled (Wardle 1990, see also Smith 1992), 
it overestimates the thickness in the present case, when the 
grains do not drift with the magnetic field. 

Numerical integration of the differential equations for- 
wards (backwards) in z from the upstream (downstream) 
state yield the shock structure, i.e B x and B y (and hence 
velocity components and fluid densities), as a function of 
z. The derivatives of B x and B y vanish far upstream and 
downstream (in the limit that \z\ — > oo), so the field com- 
ponents have to be perturbed before commencing the in- 
tegration. Because the Hall parameters are extremely large 
for the ions and electrons (~ 10 4 , and 10 6 respectively) the 
equations are stiff (this could be removed by assuming that 
the ions and electrons are well-tied to the field lines), so I 
use the integration method of Gear (1971). 



4 RESULTS 

In this section, I present solutions for C-type shock waves 
propagating into dense molecular gas. Initially, a single-size 
grain model is adopted, with grain radius 0.4 /i and abun- 
dance n 9 /nH = 2.5 x 10 -14 . The pre-shock density is set to 
uh = 10 6 cm -3 , Bq = 1 mG, and 6 = 30°, and the pre-shock 
ion density is ni/nn — 10~ 8 . 

Mildly super- Alfvenic shocks with v s /va = 1-5 are con- 
sidered in §4.1. There is a family of intermediate shocks and 
a single fast shock solution at this shock speed. Stronger 
shocks, with (v s /va = 10), for whic h only the fast shock so- 
lution exists, are considered in §4.2. For the chosen parame- 
ters, the fast shock is significantly non-coplanar. The effect 
of adopting a grain size of 0.1 fj., with a larger grain abun- 
dance, and of increasing the pre-shock density to 10 7 cm~ 3 
are also examined. 

Finally, the dynamical effects of an MRN grain-size dis- 
tribution are examined in §p^3, where solutions are presented 



for cases with and without PAHs present. 



4.1 Weak fast and intermediate shocks 

For these solutions we adopt 100 and 200 K as the effec- 
tive values of neutral and electron temperatures respectively, 
consistent with the results of Pilipp & Hartquist (1994), and 
the effective grain-neutral drift speed is set at 1 kms -1 . 
With these choices, the mean grain charge Z g — —19.15, 
and the ion, electron and grain Hall parameters in the pre- 
shock gas are 4600, -2.503 x 10 6 , and -0.1606 respectively. 

Both intermediate and fast shocks exist for these param- 
eters, with B x /Bq equal to 0.5 for the upstream state, and 
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Figure 4. The phase space for shocks with v s /va = 1.5 and 
= 30°. Each shock solution is represented by a trajectory 
(B x (z), B y (z)). The point U indicates the upstream state, and 
the points F and I indicate the downstream states for fast and 
intermediate shocks respectively. The solid curve running from 
U to F is the trajectory for the fast shock solution, and the two 
solid curves running from F to I are intermediate shocks that can 
exist in the downstream state of the fast shock. The remaining 
trajectories (dotted and dashed lines) are a one-parameter family 
of intermediate shock waves. The structure of the three interme- 
diate shock solutions indicated by dashed trajectories are plotted 
in Figures §, and |. 



1.5 and —2.0 for the downstream states of fast and interme- 
diate shocks respectively. Before considering the structure of 
individual shock waves it is instructive to examine the en- 
tire set of shock transitions in the B x -B y phase space (see 
§2.3). This is presented in Figure ^j, with points marked on 
the B x axis corresponding to the upstream state (U) and to 
the downstream states for fast (F) and intermediate shocks 
(I). Only integral curves corresponding to shock solutions, 
i.e. beginning and ending at stationary points, are plotted . 

The one-parameter family of solutions running from U 
to I are intermediate shocks, in which B± rotates by 180° 
through the shock front. A unique solution, the fast shock, 
runs from U to F close to the B x axis. In addition, two in- 
termediate shocks run from F to I with opposite senses of 
transverse field rotation. There are no other acceptable so- 
lutions, as the remaining curve terminating at F, and the 
remaining set of curves terminating at I are physically un- 
acceptable, with B diverging at z = — oo. 

The slight asymmetry between the half-planes B y > 
and By < reflects the asymmetry in the proper- 
ties of charged particles of different sign, namely that the 
negatively-charged grains are partially decoupled from the 
magnetic field, whereas the remaining charged species are 
well coupled. If all charged species were well coupled, the 
diagram would be symmetric upon reflection about the B x 
axis. 

Figure ^ and ^| show the detailed structure of two of the 
intermediate shocks whose trajectories in Fig. ^ are indi- 
cated by dashed lines. The rotation of Bx has the opposite 




12 3 4 

z (10 16 cm) 

Figure 5. The structure of the intermediate shock corresponding 
to the dashed trajectory leaving the point U in figure U towards 
the upper half-plane of phase space. The top panel shows the 
transverse velocity components of the neutrals (solid), ions and 
electrons (dashed), and grains (dotted). The z-component of ve- 
locity is plotted in the middle panel, and the transverse magnetic 
field components are plotted in the lower panel. 



sense in each shock. The shocks are similar to those pre- 
sented by Pilipp & Hartquist (1994) (see their Figs la' -lc' 
and 2a'), with B x changing sign across the shock front, and 
the {/-components of velocity and magnetic field becoming 
comparable to the ^-component within the shock front. The 
difference between the shock structures lies primarily in the 
signs of the y-components. This would be the only difference 
if all of the charged species were tied to the field lines, when 
the phase space is symmetric under reflection in the a;-axis. 
However, as the grains are poorly coupled to the magnetic 
field the grain drift through the neutrals is much smaller 
than that of the ions, electrons and magnetic field, and the 
direction of this drift shows distinct differences between the 
two shocks. For example, the grains are declerated before 
the neutrals in the middle panel of Fig. [| and lag the neu- 
trals in the middle panel of Fig. ^| Similar differences exist 
in the x and y velocity components. 

Although these two shocks are fairly typical members 
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Figure 6. As for Figure g, but for the intermediate shock whose 
dashed trajectory leaves U in Figure [l towards the lower half- 
plane of the phase space. 



Figure 7. As for Figure g, but for the intermediate shock whose 
dashed trajectory runs through the origin of the phase space of 
Figure W. Note the rarefaction precursor (z < 2 X 10 16 cm). 



of the intermediate family, shocks with trajectories pass- 
ing close to the origin of phase space cannot be charac- 
terised by rotation of B±. The structure of the interme- 
diate shock which has a trajectory passing through the ori- 
gin of phase space (see Figure |^) is plotted in Figure 0. The 
transverse field does not increase monotonically through the 
shock front, but drops initially. As the sum of the trans- 
verse magnetic field B^/8tt and the ram pressure pov s v z is 
constant through the shock front, there is a compensating 
increase in the z components of the fluid velocities, and a 
corresponding drop in density until compression begins in 
the shock (at z^2x 10 16 cm). The shock has a rarefaction 
precursor, which sucks the upstream material towards the 
shock before accelerating and compressing it. 

Rather than show the fast shock running from U to F 
and those running from F to I seperately, Figure ^ shows the 
structure of an intermediate shock whose phase-space tra- 
jectory lies close to the fast shock trajectory, almost reach- 
ing F before turning and running towards U close to the F 
to I transition with B y > 0. Because dB x /dz and dB y /dz 
become small in the neighbourhood of F, the shock struc- 
ture consists of two distinct parts, corresponding to a fast 



shock (z^6x 10 16 cm), followed by an intermediate shock 
(z i> 8 x 10 16 cm). By choosing trajectories that approach 
successively closer to F before turning towards I, the gap be- 
tween the transitions can be made arbitrarily large, becom- 
ing infinite in the limit that the trajectory actually reaches 
F. Alternatively, choosing trajectories that turn away at suc- 
cessively earlier stages causes the two structures to merge. 

The structure of the fast shock is very similar to 'stan- 
dard' solutions, with small departures of vector quantities 
out of the x-z plane. This is because the grains are not dom- 
inating the dynamics of the shock front. The z-component 
of the drift velocities is relatively small in the intermediate 
part of the shock transition because the compression across 
this part of the structure is only 4/3, compared to 3 in the 
fast part of the transition. 



4.2 Strong C-type fast shocks 

Now consider stronger shocks with v s /va = 10, correspond- 
ing to a shock speed of 18 kms -1 . As these shocks are 
characterised by larger drift speeds and higher tempera- 
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Figure 8. As for Fig. |fj, but for an intermediate shock whose 
trajectory almost reaches the point F in Fig. [j| before turning 
towards I. The transition from upstream to downstream is made 
in two distinct stages, a fast shock transition (2 < 6 X 10 16 cm) 
followed by an intermediate shock (z > 8 X 10 16 cm) (see text). 



Figure 10. As for Fig. ^| but for the fast shock of Fig. |9j. Note 
that the curves in the top panel terminating at 0.1 and are the 
x and ^-components of the fluid velocities respectively. 




5 10 15 

B x /B 

Figure 9. Integral curves in the B x -B y phase space for v s /v a = 
10, and 9 = 30°. The solid curve running from U to F is the 
fast shock trajectory. The neighbouring trajectories (dashed) all 
correspond to unphysical shock solutions. 



tures, I adopt 1000 K, 2000 K, and 10 kms" 1 as the ef- 
fective values of the neutral and electron temperatures and 
the grain drift speed. This implies a mean grain charge 
Z g = —191.5, and ion, electron, and grain Hall parameters 
of 4609, -7.915 x 10 5 , and -0.2420 respectively. 

The shock speed is too high for intermediate shocks 
to exist as v s /va > V2cos0 « 2.45. The fast shock running 
from U to F remains viable. The structure of the phase space 
around the fast shock trajectory is shown in Figure ^. Each 
of the neighbouring trajectories peels off and heads towards 
I, eventually becoming unphysical. Thus an attempted inte- 
gration of the fast shock solution starting at U will stumble 
onto a trajectory headed for I, and yield only unphysical 
solutions, as found by Pilipp & Hartquist (1994). However, 
integration in the reverse direction, from F to U, will be 
successful. 

The resulting fast shock solution, shown in Figure |l^, is 
similar to a 'standard' C-type shock, apart from departures 
of the magnetic field and fluid velocities out of the x-z plane 
at the level of 50 per cent. The grains lag the ions in the x 
and z directions, where they are caught between the mag- 
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Figure 12. The heating rate per unit mass within the shock 
fronts of figures lu^ (peaks at z ft: 5 X 10 15 cm) and [ll] (peaks 
at z £« 6.5 X lO^cm). For each shock the heating associated 
with grain-neutral, ion-neutral, and electron-neutral collisions are 
plotted as solid, dashed, and dotted curves respectively. 
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Figure 11. As for Fig. |l^, but with dB y /dz set to zero. 



netic field and neutral fluid; and they lead in the y-direction 
in which they are responsible for the excursion. 

To evaluate the effect of the out-of-plane forces on the 
shock structure, the result of supressing B y is plotted in Fig- 
ure [ll]. The y-component of the ion and electron velocities 
are very small, as these species are tied to the field lines. 
The grains, however still have a substantial drift in this di- 
rection, but this component of the drag force has been ne- 
glected. The shock is substantially thicker (5 x 10 15 cm) than 
the full solution shown in Figure [Hi] (3 x 10 15 cm), and drift 
velocities are much higher in the full solution as the kinetic 
energy flux that must be dissipated into thermal energy is 
the same in both cases. The temperature within the shock 
of Fig. |l(^ will be significantly higher. Although the temper- 
ature cannot be calculated in our prescription, it is useful to 
examine the heating rate per unit mass associated with the 
collisions between neutrals and charged species, 



— =2^ l liPi\vi -v\ 



(81) 



(Draine 1986; Chernoff 1987), which would largely be bal- 
anced by the local radiative cooling rate per unit mass. JjPjo 



is 1.485 x 10" 



8.649 x 10" 



1.354 x 10" 



for the 



ions, electrons and grains respectively. Thus for a given drift 
speed, the grains contribute a drag force and collisional heat- 
ing rate ten times that of the ions, and the heating by elec- 
tron collisions is negligible. The actual heating rate for the 
different species within the two shocks is plotted in Figure 
|l2} As expected, the heating is more severe in the full solu- 
tion. However, note that most of the difference is due to the 
grain collisions which dominate the ion collisions in the full 
solution. It is clear, therefore, that the artificially coplanar 
solutions that have been used in detailed C-shock models 
(Draine et al. 1983; Kaufman & Neufeld 1996a, b) may sub- 
stantially understimate the the peak shock temperature. 

Now consider the effect of raising the pre-shock density 
to 10 7 cm -3 . This increases the importance of the grains, as 
the pre-shock ion fraction is decreased by a factor of three, 
and the grains are further decoupled from the magnetic field. 
The pre-shock ion, electron, and grain Hall parameters in 
this case are 460.9, -7.915 x 10 4 , and -0.02420 respec- 
tively. The solution is plotted in Figure |l3[ As expected, 
the y-components of the fluid velocities and magnetic field 
are increased, and the grain drift through the neutrals is 
substantially decreased. There is also the the beginnings of 
oscillatory behaviour in the shock precursor, corresponding 
to the spiralling of the shock trajectory out of the point U 
in the shock plane (see Fig. |l6| ). 

Finally, consider the effect of decreasing the grain size 
from 0.4 to 0.1/i, which increases the grain abundance by 
a factor of 64 to n g /nu = 1.6 x 10 -12 . As the mean grain 
charge is reduced by a factor of 4, to Z g = —47.87, this 
increases the charge resident on grains by a factor of 16, 
and grains now dominate the current within the shock front. 
The shock solution is plotted in Fig. jw| Although the ion- 
neutral streaming velocities are a significant fraction of the 
shock speed in the precursor, most of the energy is dissipated 
by grain-neutral collisions within the extremely thin main 
part of the shock front (z = 0.14-0.18 x 10 15 cm) - ions 
would have to stream through the neutrals over a length 
scale of order 5 x 10 15 cm to dissipate the incoming kinetic 
energy. Besides producing significant partial rotation of the 
field and fluid velocities within the shock front, the precursor 
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of the shock exhibits an oscillatory behaviour because the 
propagation of wave modes in the pre-shock gas at these 
length scales is substantially modified. This is manifested in 
the linear analysis near the stationary points, which shows 
that the upstream state becomes a strongly spiral node. The 
ion and electron velocities approach zero and their densities 
rapidly increase at z ~ 0.15 x 10 15 cm, at which point large 
currents prevent further movement to the left in the B x -B y 
phase space (see Fig. |l6| ) and the B x and v x profiles are 
clipped. The field-line structure of the shock is sketched in 
Fig. @ 

The trajectories of the shock solutions are plotted in 
Figure |l^. As expected, the shocks are no longer close to 
being coplanar, and the upstream state becomes a spiral 
node, with the magnetic field oscillating ahead of the shock 
solution. The long-dashed curve shows the locus where the 
the ion and electron densities diverge, i.e. where Vi Z = v ez = 
(see §3.4). The trajectory for the shock propagating into 



the 10 7 cm a gas with 0.1 /i grains approaches this locus, 
but is forced away by the associated large current density. 



shock front 



shock velocity 




x y 



X Z X z 

face-on side 



Figure 15. The magnetic field structure in the shock solution 
plotted in Fig. [l4| 
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Figure 16. The phase portraits for the fast shock solutions of 
Figures [l(^ (solid), ^ (short-dashed) and ^ (dotted). The long- 
dashed curve indicates where the ion and electron densities di- 
verge (see text). 



4.3 The effect of the grain-size distribution 

In this subsection, I consider the effect of an MRN distri- 
bution on the shock dynamics, both with and without an 
additional population of very small grains. Ionization mod- 
els of the ambient molecular gas show that at gas densities 
of order 10 6 cm -3 , not all grains are charged, but the no. 
density of those that are is order 10 _11 nH- Thus, the net 
charge on the grains, n g Z g « 10 _11 nH is of the same order 
as that for the single-grain-size models presented in the pre- 



vious subsection, for which n g /nn ~ 10~ and Z g « —120. 
The shock models are therefore rather similar. 

When PAHs are present, they are the dominant charged 
species. The effects of interest are greatest at densities 
around 10 6 cm -3 , as at higher densities the abundances of 
positively and negatively charged grains become very sim- 
ilar, and |<7i|/cr2 becomes small. For tih = 10 6 cm -3 , the 
fractional abundances of positively charged PAHs plus ions 
is 3 x 10~ 9 , with a similar abundance of negatively charged 
PAHs. Singly-charged grains and electrons both have abun- 
dances of approximately 2 x 10 -11 (Nishi et al. 1991; Kauf- 
man & Neufeld 1996a). Thus I adopt pre-shock densities 
and magnetic fields tih = 10 6 cm -3 and Bo = 1 mG and a 
four fluid model consisting of neutrals, positively and nega- 
tively charged ion-like particles each with n/nn = 2 x 10 , 
negatively charged grains with N/nu = 2 x 10 -11 . The 
two species of ion-like particles are assumed to be iden- 
tical apart from the sign of their charge, with < crv > = 
1.6 x 10~ 9 cm~ 3 s _1 , so their pre-shock Hall parameters are 
±4609. The grains are assumed to satisfy the power-law size 
distribution ( p6] ) between a\ = 50A and = 2500A. Given 
a choice of the pre-shock Hall parameter of the smallest 
grains, Pi, the effective values of n g and j3 g can be calcu- 
lated. Rather than show the detailed shock structure, which 
is similar to that presented previously, I show in Figure [l?] 
the phase-space trajectories for different choices of /3i. As /3i 
decreases, f3 g decreases, and n g increases, so the shocks be- 
come increasingly non-coplanar, with the characteristic spi- 
ral trajectory for \/3 g \ ^ 0.3. This is to be expected from the 



Figure 17. Integral curves in the B x -B y phase space for shocks 
with njj = 10 6 cm" 3 , v s /va = 10, and 9 = 30°, for the case when 
grain sizes follow an MRN distribution and there is an additional 
population of PAHs. The different trajectories are labelled by the 
value of the pre-shock Hall parameter for the smallest grains, with 
radii of 50 A (see text). 



non-coplanarity criterion ([t5|). These results indicate that 
significant non-coplanarity occurs for B\ 1, that is, in re- 
gions of unusually weak field (see eq |p8[ ). 

In the absence of PAHs, the dominant charged fluids 
for nn 5i 10 9 cm -3 are ions, electrons and grains carry- 
ing a single negative charge (Nishi et al. 1991). At den- 
sities of 10 6 , 10 7 , 10 8 cm , I adopt fractional abundances 
n,/n H = 10~ 9 , 3 x 10~ 10 , 1.03 x 10~ 10 and jV/n H = 3 x 
10~ 10 , 2 x 10 -10 , 1 x 10 -10 respectively. At the two lower 
densities, the density of electrons is determined by charge 
neutrality once the charged grain abundance is corrected 
for the size distribution. At 10 8 cm -3 , the electron density 
is relatively small and would be strongly affected by this pro- 
cedure, so the no. density for charged grains is not corrected 
for the size distribution. In any case, the corrections would 
amount to changes of order 5 per cent. Figure [l8| presents 
the trajectories for shocks with v 3 /va ~ 10 at these three 
densities for standard and lower values of the magnetic field. 
This illustrates how the shock trajectory depends primarily 
on the ratio B/nu that determines the Hall parameters, with 
the particle abundances and field strength determining the 
shock thickness. 



5 DISCUSSION 

The results presented in §^ confirm that both intermedi- 
ate and fast C-type shock solutions exist. A one-parameter 
family of intermediate shocks exists for each shock speed 
between the intermediate speed va cos 9 and V2va cot 9; at 
higher speeds the downstream states become unphysical. A 
unique fast shock exists for each shock speed in excess of the 
fast speed va- The relationship between these two classes 
of solution is summarised by the trajectories in the B x -B y 
phase-space diagram. 

The fast shock trajectory runs along the B x axis if the 
charged species are tied to the magnetic field lines, or if 
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Figure 18. Integral curves in the B x -B y phase space for 
shocks with v s /va = 10, and 9 = 30°, for the case when 
grain sizes follow an MRN distribution and PAHs are absent. 
The different trajectories correspond to different choices of pre- 
shock magnetic field and gas density — solid:l mG, 10 6 cm -3 ; 
short-dashed:0.3mG, 10 6 cm -3 or 3mG, 10 7 cm , long- dashed: 
lmG, 10 7 cm" 3 or 10 mG, 10 8 cm- 3 ; dotted: 3mG, 10 s cm" 3 . 



positive and negatively charged species have identical abun- 
dances and Hall parameters. In this case, integration of the 
fast shock solution away from the upstream state is straight- 
forward. The direction of the trajectory leaving the up- 
stream stationary point (i.e. along the B x axis) is known, 
and the velocity and magnetic field vectors in the fast shock 
solution are confined to the x-z plane, so the y components 
can be discarded. Integration can be started by stepping 
off U towards F, and then integrating the trajectory (e.g. 
Draine 1980; Wardle & Draine 1987). 

The symmetry of phase space under reflection in the B x 
axis is destroyed by grains, which being partially decoupled 
from the magnetic field and preferentially carrying a nega- 
tive charge, introduce a handedness into the fluid. The fast 
shock trajectory no longer runs along the B x axis, and this 
distortion provides an opportunity for numerical mixing of 
the fast solution with the neighbouring divergent interme- 
diate shock trajectories during integration. Any attempted 
integration from U to F is doomed. To begin with, there is 
no a-priori means of determining the direction for the initial 
step off the point U. In principle this problem could be cir- 
cumvented by using a shooting method to adjust the starting 
direction until the integration successfully reached F. How- 
ever, because of finite numerical precision, any attempted 
integration into F along the fast trajectory will step onto a 
neighbouring intermediate shock trajectory and eventually 
be forced away from the fast trajectory and towards I (c.f. 
integration towards sonic points) . This behaviour prevented 
Pilipp & Hartquist (1994), who started their integrations at 
U and treated the direction of the first step away from U as 
a free parameter, from finding the fast shock solutions. At 
low shock speeds, they discovered the family of intermediate 
shocks but could not find the fast shock, noting instead that 
the sense of rotation of B± within the intermediate shocks 
changed at a critical direction of the initial step away from 
U. At higher speeds, for which the intermediate shock solu- 



tions become unphysical, Pilipp & Hartquist were unable to 
find any acceptable solutions. 

A realistic treatment of the ionization balance is re- 
quired for a quantitative assessment of the role of grains in 
C-type shocks. The ionization balance is particularly impor- 
tant when v z for any charged species approaches zero and 
the number density of that species diverges in the constant- 
flux approximation (eq. ^|) used here. This does not affect 
the divergence of trajectories from the fast trajectory un- 
less the fast trajectory itself approaches this limit (e.g. the 
trajectory shown in Fig. [lfi] for the shock displayed in Fig. 

Some conclusions are robust to changes in the micro- 
physics. If the conductivity of the gas can be regarded as 
some unspecified function of local physical variables at any 
point within the shock front, so that shock solutions can still 
be represented as trajectories in the B y -B x phase-space, 
the topology of the phase space is determined by the na- 

to be 



3.2 



ture of the fixed points. These were shown in 
a source, sink and saddle independent of the conductivity 
tensor. Thus the integral curves always have the topology 
plotted in Figure ^, apart from smooth distortions (as appar- 
ent in Fig ^) because of changes to the conductivity tensor. 
The intermediate-shock curves always run from U or D to I, 
and thus Vi Z = v ez = at some point if the shock speed is 
high enough for intermediate shocks to be excluded. The fast 
shock trajectory could in principle cross the Vi Z — v ez = 
locus if the ionization balance is able to distort the trajec- 
tory sufficiently, although it is apparently unable to do so in 
the fixed-flux approximation used here (see Fig. p"s| ) . In any 
case, when the fixed-flux condition is dropped, crossing this 
locus does not imply infinite charged-particle abundances 
as the increase is limited by recombinations. Intermediate 
shocks, however, will still be excluded as their downstream 
conditions are unphysical. 

A more realistic treatment of the physical processes oc- 
curring within the shock front is essential as the grain prop- 
erties within the shock and the shock structure itself are 
intimately related. In particular, the neutral and electron 
temperatures are required to accurately calculate the Hall 
parameters for a given grain size. In addition, the grain-size 
distribution varies within the shock front as different sized 
grains will generally have different drift speeds, the smallest 
grains being compressed with the magnetic field, the largest 
grains with the neutrals. Incorporation of these effects in- 
creases the dimensionality of the phase space by introduc- 
ing differential equations for the temperatures and relative 
abundances of different species. In principle this could alter 
the topology of the phase space, although this is unlikely to 
be the case in practice. There are, however, practical prob- 
lems in finding the fast shock structure by integrating back- 
wards from the downstream state if chemistry is coupled to 
the hydrodynamics (through e.g. changing the abundances 
of coolants or the degree of ionization in the shock), since 
one does not then know the post-shock conditions a-priori 
(c.f. the discussion in Roberge & Draine (1990) regarding 
C*- and J-shocks), and a time-consuming shooting method 
must be employed. 

Although the grain drag significantly modifies the shock 
structure, it is unclear whether this will produce gross dif- 
ferences in the total line emission from the shock front. The 
essential feature of C-shock structure is that it maintains 
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a layer of hot molecular gas. For a given pre-shock den- 
sity and shock speed, the same power is dissipated per unit 
area independently of the detailed structure of the shock 
front. The gas generally will be heated up to temperatures 
between 1000 and 2000 K at which point the local heat- 
ing and cooling balance. Molecules will remain intact and 
the molecular line emission will not be greatly affected. The 
insensitivity of line ratios to details of shock structure is il- 
lustrated by recent numerical work (Stone 1997; Neufeld & 
Stone 1997; MacLow & Smith 1997), which follows the de- 
velopment and saturation of the instability in C-type shocks 
(Wardle 1990). Even though the neutral gas is collected into 
dense fingers within the shock front, the line ratios are gen- 
erally unchanged by more than a factor of two (Neufeld & 
Stone 1997), although diagnostics can be found (MacLow & 
Smith 1997). Nevertheless, one might expect that the grains 
will affect the stability of C-shocks and also the speed at 
which shocks become J-type. 

Finally, it is worth pointing out that the theoretical 
study of C-shocks in molecular clouds will contribute to the 
general theory of MHD shock waves, and of intermediate 
shocks in particular (see the review by Wu 1995). Until re- 
cently intermediate shocks in MHD were thought to be un- 
physical because they appeared to lack sufficient freedom 
to adjust to slight perturbations in the upstream or down- 
stream flow (see, e.g. Kantrowitz & Petschek 1965). The is- 
sue was reopened when intermediate shocks emerged as sta- 
ble structures in numerical simulations (Wu 1988a; Steinolf- 
son & Hundhausen 1990a, b). It has since been realised that 
the arguments against the existence of intermediate shocks 
are based on ideal MHD, which breaks down completely 
within the shock front. Numerical studies of time-dependent 
intermediate shocks in resistive MHD demonstrated that the 
shocks could be formed 'naturally' and were stable (Wu 
1990). An examination of the structure of resistive, inter- 
mediate shocks in the weak shock limit (Kennel et al. 1989; 
Wu & Kennel 1992), shows that, as found for C-type shock 
waves, for given external parameters (i.e. pre-shock condi- 
tions and shock speed), a one-parameter family of interme- 
diate shock structures exist connecting the jump conditions. 
This 'internal' parameter provides the additional freedom for 
the shock to adjust its structure to external perturbations, 
a behaviour outside of the scope of ideal MHD in which the 
shock transition is represented as a discontinuity between 
the upstream and downstream states. As these conclusions 
depend on shock structure, in principle they depend on the 
model used for the magnetised medium. Resistive (Wu 1990) 
and hybrid (Wu & Hada 1991) models have been investi- 
gated so far. With the advent of multidimensional ambipolar 
diffusion codes, it has become possible to study these issues 
in weakly ionized media (see Smith & MacLow 1997). 



6 SUMMARY 

In this paper I constructed steady models of oblique C-type 
shocks in which the magnetic field within the shock front is 
not artificially confined to the plane containing the upstream 
and downstream magnetic field and the shock normal. Four 
fluids - neutrals, ions, electrons or PAHs, and negatively 
charged grains - were considered. The thermal pressure of 
the fluids and the inertia of the charged components were 



neglected. The effects of chemistry and changes in fractional 
ionization, and the charge residing on grains were ignored. 
The rate coefficients for grain-neutral and electron-neutral 
elastic scattering (which depend in principle on the neutral 
temperature and the grain-neutral drift speed, and on the 
electron temperature, respectively) were assumed constant 
with values appropriate to conditions within the shock front. 

These assumptions reduce the number of differential 
equations describing the shock structure to a pair for the 
components B x , B y of the magnetic field that are perpen- 
dicular to the shock normal. Shock solutions can there- 
fore be conveniently represented by trajectories in a two- 
dimensional phase space. 

Models were presented for weak shocks (v s /va = 1-5), 
and for strong shocks (v s /va = 10). Variations in the grain 
size from 0.1/x to OAfj,, and the effect of an MRN grain-size 
distribution were considered. 

The results are summarized as follows: 

(i) The cold MHD jump conditions permit three station- 
ary points in the phase space, corresponding to the upstream 
state and downstream states of the fast (for v s > va) and 
intermediate shocks (va cos 9 < v s < VACot9). Valid shock 
solutions have phase-space trajectories linking the upstream 
state to one of the downstream states. A linear analysis, 
valid for any number of charged species, shows that the up- 
stream state is a source, the fast downstream state is a sad- 
dle, and the intermediate downstream state is a sink. 

(ii) A one-parameter family of intermediate shocks exists 
for each shock speed in the range vacosO < v s < vacoXO 
where 9 is the angle between the shock normal and the pre- 
shock magnetic field. These solutions correspond to those 
found by Pilipp & Hartquist (1994). The family contains 
members with either sense of rotation of B± through the 
shock front, members with a rarefaction precursor in which 
B± becomes small within the front before compression be- 
gins, and members that correspond to a fast shock followed 
downstream at a finite distance by an intermediate shock. 

(iii) A unique fast shock exists for each shock speed 
v s > va, where va is the Alfven speed in the pre-shock gas. 
Pilipp & Hartquist (1994) were unable to find these solutions 
because integration of the equations for shock structure from 
the upstream state to the downstream state is unstable. In- 
stead, integration must begin at the downstream state and 
run backwards through the shock front. 

(iv) When all charged particles are well couwell coupled 
to the magnetic field, or when there is a symmetry between 
the poorly-coupled particles of either sign, the phase space 
is symmetric about the B x axis, and the fast shock is copla- 
nar, its trajectory running along the B x axis. In conditions 
typical of molecular clouds, negatively-charged grains con- 
tribute significantly to the drag and are loosely coupled to 
the field. The net asymmetry in the coupling of negatively 
and positively charged particles to the magnetic field lines 
imposes a handedness on the shock structure and the phase 
space trajectories lose their symmetry on reflection in the 
B x axis. In particular, the fast shock trajectory no longer 
runs along the axis - the magnetic field and fluid velocities 
no longer lie in the x-z plane (containing the pre-shock field 
and shock normal) but have significant y components within 
the shock front. 

(v) For typical conditions in molecular clouds, grains 
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dominate the frictional force and the magnetic field is well- 
tied to ions, electrons and PAHs. Under these conditions 
for given values of pre-shock magnetic field and gas den- 
sity, the pre-shock grain Hall parameter determines the 
shock structure apart from an overall length scale for the 
shock thickness which is determined by the grain abundance. 
Thus under the approximations adopted in this paper, the 
shock structures are relatively insensitive to the details of 
the grain-size distribution and the exact composition of the 
charged species in the gas. 

(vi) An MRN grain-size distribution can be approxi- 
mately incorporated by calculating an effective abundance 
and an effective grain Hall parameter for single-size grains. 
The ratio of the effective abundance to the abundance of 
charged MRN grains, and the effective Hall parameter are 
determined by the Hall parameter of the smallest grains in 
the MRN distribution. 

(vii) The degree of noncoplanarity is determined by the 
grain Hall parameter (3 g , with significant effects being found 
for \p g \ SS 1. When |/3 S | ^ 0.3 the upstream state becomes a 
spiral node and the shocks exhibit a precursor in which B± 
may make several rotations. 

(viii) Ions and electrons become passive in the sense 
that they may stream through the neutral gas with large 
transverse drift velocities (of order half the shock speed 
or more) without generating significant dissipation, as the 
shock thickness is determined by the collisions of neutrals 
with charged grains. The primary role of these charged 
species is to guarantee that the electric field in the shock 
front is very nearly orthogonal to the magnetic field. 

(ix) Supressing the out-of-plane components of the drift 
velocities and magnetic field may significantly reduce the 
current and leads to a thicker shock structure (by, e.g. a 
factor of two) and a consequent decrease in the heating rate 
within the shock front. Models supressing the out-of-plane 
components of the magnetic field and velocity therefore un- 
derestimate the temperatures within the shock front, and 
significantly underestimate the magnitude of the grain drift 
speed through the neutrals. 
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